Collective Thermotaxis of Thermally Active Colloids 
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Colloids with patchy metal coating under laser irradiation could act as local sources of heat due 
to the absorption of light. While for asymmetric colloids this could induce self-propulsion, it also 
leads to the generation of a slowly decaying temperature profile that other colloids could interact 
with. The collective behavior of a dilute solution of such thermally active particles is studied using 
a stochastic formulation. It is found that when the Soret coefficient is positive, the system could be 
described in stationary-state by the nonlinear Poisson-Boltzmann equation and could adopt density 
profiles with significant depletion in the middle region when confined. For colloids with negative 
Soret coefficient, the system can be described as a dissipative equivalent of a gravitational system. 
It is shown that in this case the thermally active colloidal solution could undergo an instability at 
a critical laser intensity, which has similarities to supernova explosion. 

PACS numbers: 82.70.Dd,05.70.Ln,47.57.-s,47.70.-n 



The motion of colloidal particles in a solution in the 
presence of an externally applied temperature gradient, 
which is known as thermophoresis or the Soret effect [ij , 
has been studied since the 19"^ century and observed 
in a variety of systems [H. Although its existence can 
be well formulated in non-equilibrium thermodynamics 
within linear response theory [3| , many aspects of the mi- 
croscopic nature of the phenomenon has remained a sub- 
ject of investigation to date The effect has also been 
shown to provide a powerful tool for manipulating macro- 
molecules and colloids Q. Since in the phoretic trans- 
port mechanisms the colloids experience no net force, 
it is possible to take advantage of them to design self- 
propelled particles by incorporating a built-in mechanism 
that provides asymmetric sources that could generate 
and maintain the necessary gradient across them needed 
for propulsion [1,01- Recently, Jiang et al. Q have shown 
that silica beads half-coated with gold when irradiated 
with a defocused laser beam exhibit such a propulsion, as 
the gold caps act as heat sources when they absorb light. 
Moreover, even without the self-propulsion, laser-heated 
gold-coated colloids have been shown to undergo sub- 
stantially enhanced Brownian diffusion, which is related 
to the modification of the temperature in the medium 
and the resulting changes in the viscosity Q. Since such 
thermally active colloids would create temperature pro- 
files around them that decay as l/r, in addition to caus- 
ing them to self-propel thermophoresis could provide a 
mechanism for them to interact with one another in a 
solution. The long-ranged nature of the inter-colloidal 
thermophoretic interaction could lead to interesting col- 
lective behaviors. 

Here we construct a stochastic formulation to describe 
the collective behavior of a number of thermally active 
colloids. At the long time and large length scale limit 
and for dilute solutions, the formulation simplifies to a 
set of two nonlinear coupled differential equations for 
the density and temperature profiles in the medium. In 



stationary-state, we provide a number of examples for 
which the equations could be solved exactly. They show 
a depletion effect for the case of positive Soret coefficient, 
and an instability at a finite laser intensity for negative 
Soret coefficient. 

We consider N colloidal particles of radius R that are 
half-coated with a metal that absorbs the laser light with 
an efficiency e, thus creating a local source of heat of mag- 
nitude el, where I is the intensity of the laser (see Fig. [IJ . 
We assume that the laser intensity is uniform through- 
out the space, and thus ignore any optical confinement 
effect. In an externally generated temperature gradient, 
the colloids move with a drift velocity v = — £>tVT, 
where Dt is the thermodiffusion coefficient. The asym- 
metric heat generation around each colloid provides a 
mechanism to create and maintain a local temperature 
gradient that leads to propulsion via self-thermophoresis 
[7|. For a Janus-sphere colloid, the propulsion veloc- 
ity can be calculated as wq = elDx/idfi), where k is 
the thermal conductivity of the medium 0, H, • The 
stochastic motion of the i-th colloid is described by its 
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FIG. 1: (color online.) The metal-coated Janus-spheres under 
irradiation could self-propel in the directions shown by the 
(green) arrows and interact with one-another via the long- 
ranged temperature profiles they generate. The interactions 
are mutually repulsive when St > and attractive when 
St < 0. 
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instantaneous position ri{t) and orientation ni{t) that 
is a unit vector. They satisfy the Langevin equations 
^ = Von, - DrVTiYi) + and ^ n, x n„ in 
which and rj^ are Gaussian-distributed noise terms. 
Since each colloid generates heat with an overall (surface- 
average) rate of \el, the temperature profile at the loca- 
tion of each colloid is affected by the heat generated by 
all the other colloids. Since the temperature profile equi- 
librates considerably faster than the colloids, we have 

r(r,t) = To + |r-r,(t)h *° ^^"^ 1°"^*^^* o^^*^^ 

the multipole expansion of the heat source distribution. 
While this approximation should be valid for sufficiently 
dilute colloidal solutions, the contributions from higher 
multipoles could readily be added to the above temper- 
ature profile. To probe the colloidal activity at shorter 
times than the rotational diffusion time, we need to in- 
corporate the time dependence of heat diffusion, which 
could lead to anomalous dynamics of the colloid jTlj. 

The Fokker-Planck equation for the probability dis- 
tribution P(r,n,<) = ^E^=i'5(r-r.(t))<5(n-n,(<))^, 
can be constructed from the Langevin equations as 

dtV+V- \vQnT~DT (VT) V-DVv] -DrU^V = 0, (1) 

where 7?, = n x (9„. In Eq. ([l}, D and are the 
translational and rotational diffusion coefficients, respec- 
tively, and represent the corresponding widths of the 
Gaussian probability distributions for the noise terms 
and rii in the Langevin equations. In a medium with 
uniform temperature T, we have D = k^T / {Q-ki]R) and 
Dr = k-gT / {^nrjR^) , where -q is the viscosity of water. 
Equation ([T|) should be complemented with the heat dif- 
fusion equation 



/ P(r,n), 

J n 



(2) 



which describes how the temperature profile is affected 
by the spatial distribution of the colloids due to their role 
as motile heat sources. Equations ^ and ([2|) should be 
self-consistently solved to obtained the probability distri- 
bution of the colloids as well as the temperature profile 
in the medium. 

Let us define the density p(r) — J^V{r,n), the polar- 
ization field p(r) = n ^(r,!!), and the nematic order 
parameter Q(r) = [nn — ^l] V{r,n). Performing 
on Eq. (HI), we can obtain an equation for the density 



as dtp + wqV ■ p — V 



0, which 



Dt (VT) p + D\7p 
is incomplete since it has a source term in the form of 
— iiqV • p, which is present due to the self-propulsion of 
the colloids. Performing J^nx Eq. ([T]), we can obtain 
an equation for the polarization field as 



Vn 

dtp + 2DrP + -fVp + voV ■ Q(r) 
-V • \Dt (VT) p -I- DVp\ = 0, 



(3) 



where Ti?n = ~2n is used. Equation ^ is also incom- 
plete as it depends on Q, and this hierarchy will continue 
for higher order cumulants. 

To make further progress, we can seek to truncate the 
hierarchy and simplify Eq. ([3]) in some approximation. 
At time scales much longer than l/D,., the time deriva- 
tive term in Eq. ([3]) is considerably smaller than 2DrP 
and can thus be ignored. For sufficiently dilute solutions, 
i.e. when pE? <^ 1, and in the absence of any external 
means that could induce polarization, such as external 
magnetic field for particles with a magnetic dipole mo- 



ment or gravity 
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we can ignore the V • Q(r) term. 



Since the time scale is much longer 1/Dr, any transient 
or initial ordering would also have decayed. Finally, at 
length scales much larger than R, the gradient terms in 
Eq. ([3]) can be neglected, and the equation can be ap- 
proximated as 



_vo_ 

6Dr 



Vp. 



(4) 



Using Eq. (g]), we can estimate that Q ~ \7p\7p, ignor- 
ing which is consistent with our assumption of a dilute 
solution. Putting Eq. back in the density equation, 
we find 



dtp-V- D^sVp + Dt{VT)p =0, 



(5) 



where D^s = D+Vq / {6Dr) is the enhanced effective diffu- 
sion coefficient for the self-propelled active colloid 0, [Tlj , 
which could also be rewritten as T'off = T* [l -I- fPe^] (fo r 
a sphere) in terms of the Peclet number Pe = vqR/D [12l [. 
Equation ^ should be solved in conjunction with the 
heat diffusion equation, which reads — V^T = 22LilB—p^ 

In stationary state, Eq. is satisfied if DcsV p + 
Dt (VT) p — 0, which can be written as Vlnp = 
— -^^VT. If we ignore the temperature dependence in 



Z^cff! and use the Soret coefficient St 
be integrated to yield 



Dt/D, this can 



p{r) = po exp < - 



ST[T{r)^To] 
(1 + fPe^) 



(6) 



Putting the above equation back in the heat diffusion 
equation yields a single nonlinear equation for the tem- 
perature profile as 



V^T: 



27TeIR^po J ST[T{r)-To] 
exp < 



(l + fPe^ 



(7) 



Equation ([7]), which is reminiscent of the Poisson- 
Boltzmann equation for electrolytes (see below), could 
be solved for the temperature profile, which then yields 
the stationary-state density profile of the colloids via Eq. 
© . Note that the Soret coefficient could be both positive 
and negative. 
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We can define an appropriate dimensionless temper- 



ature as 
scale 



\St\IT{v)~To] 



and a characteristic length 



£IR^\St\ 



2k 1 



(8) 



which is reminiscent of the Bjerrum length in the elec- 
trostatic analogy. We have 



(9) 



where = 4:Tt£pq, and the sign choice is — sgn(S'T)- 
Equation ^ is subject to the constraint N ~ pQ e^* . 

We now discuss a number of interesting exact solutions 
for the stationary-state described by Eq. We con- 

sider colloidal solutions that are confined in some region 
in space and examine the effect of the dimensionality of 
the confinement geometry as well as the nature and the 
strength of the thermotactic coupling. The confinement 
could in practice come from the trapping effect of nonuni- 
form laser beams, which could provide a very powerful 
tuning parameter. For simplicity, we model the confine- 
ment by introducing sharp boundaries such as confining 
walls. 

When > 0, the electrostatic analogy is complete 
as the colloidal particles mutually repel one another, and 
the overall heat fiux coming out of the solution through 
the boundaries of the confining "cage" is reminiscent of 
the electric field fiux lines, which could be thought of as 
an outer shell of opposite charges maintaining neutrality 
and stability. We denote the St > case as thermo- 
repulsive. This problem can be solved exactly for ID 
and 2D confinements, and numerically for the 3D case 
15|. When the colloidal solution is confined between 
two plates of lateral size L and distance 2ft,, the density 
profile of the colloids is given as 



p{x) 



Pa 



1 



fe2 



COS- 



(10) 



(wall), and k satisfies f 



where pq is the concentration at the edge of confinement 

tan(^) = The 
profile of Eq. pi))) describes an accumulation of the col- 
loids near the confining boundary (that is reminiscent 
to counterion condensation [l^ ) and a corresponding 
depiction of the central region of the system. In the 
strong coupling limit when Nih/L'^ :s> I, we can find 
an approximate solution to the transcendental equation 



1 



In this limit, the ratio he- 



as fc/i ~ ^ 

tween the density of the colloids in the middle and at the 
edge can be found as ^ ~ i {Nih/ L'^)~'^ , which shows 
a significant depletion effect. Note that the depletion be- 
comes stronger as h is increased, when other parameters 
are kept fixed. The length 1/[2tt£{N/L'^)] is equivalent 
to the Gouy-Chapman length in the electrostatic analogy 
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FIG. 2; (color online.) The density of colloids in the middle 
of the confined space relative to the density at the edge as 
a function of the dimensionless thermophoretic coupling con- 
stant for the thermo-attractive {St < 0) case. The critical 
density at the onset of instability is (/5m/po)c ~ 3.29, which 
occurs at {AvlhN / L^)^ = 2. 



14| . For a colloidal solution trapped in a cylindrical cage 
of length L and width 2h, the density profile reads 



p{r) 



Po 



[1 



' Ni 
V L 



)Y[l-lk^r^y 



(11) 



where kh 



2+(JV£/L) • I'l'^ strong coupling limit in 
this geometry corresponds to Ni/L ^ 1. in which case 
we have ^ ~ 4 (Ni/L)^^. Note that the magnitude 
of depletion is independent of the confinement size in 
this geometry. The ratio iN/L is analogous to the so- 
called Manning- Oosawa parameter for highly charged 
rodlike polyelectrolytes [ij]. A similar profile can be 
found when the colloidal solution is confined to a spher- 
ical cage of diameter 2/i, where in the strong coupling 
limit that corresponds to N£/h 3> 1 in this case, we have 
^ ~ 21.4 {Ni/h)""^. Here, the depletion is inversely re- 
lated to the size of the cage, namely it decreases for larger 
confinement sizes. 

When 5"^ < 0, the colloids attract each other and the 
problem is analogous to a gravitational system. We thus 
denote this case as thermo-attractive. Let us go back to 
the ID confinement geometry, where the relevant ther- 
mophoretic coupling constant is Nlh/L"^ as discussed 
above. In this case, Eq. ([9]) (with the positive sign 
choice) can be integrated in closed form and the den- 
sity profile can be calculated. The stationary-state den- 
sity profile (not presented here for brevity) shows that 
the particles will accumulate towards the center of the 
confined area. Figure [2] shows that the ratio between 
the density in the middle and at the edge of the con- 
finement region increases as the thermophoretic coupling 
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constant increases, up to a critical point beyond which 
a stable (stationary-state) solution no longer exists. The 
onset of instability occurs at {'iTrehN/L'^)c = 2, at which 
{Pm/ Po)c = 3.29. Similar instabilities exist in the 2D and 
3D confinement cases [I^, where N£/L and Ni/h play 
the role of the thermophoretic coupling constant, respec- 
tively. 

The instability occurs because the particles that act 
as heat sources attract each other and could result in 
a suspension that becomes increasingly denser and hot- 
ter. In this case, the heat flux at the outer boundary of 
the system cannot balance the heat generated inside the 
confined region, which leads to an uncontrolled buildup 
of thermal energy. In molecular systems, an equation of 
the form of Eq. ([7]) (with St < 0) is used to describe 
exothermic combustion reactions that could lead to ther- 
mal explosion [l^. It is not possible to predict what 
happens in our colloidal system above the onset of this 
instability using the present formulation, as the approxi- 
mations used in its derivation are no longer valid. Within 
our handwaving analogy to a gravitational system, this 
explosion would have similarities to a type I supernova 
for a white dwarf, for which accreted material from the 
surroundings accelerates exothermic nuclear reaction to 
the point that the system becomes unstable 17|. This 



analogy is very rough, however, as the colloidal system 
operates in the dissipativc regime as opposed to the iner- 
tial and conserved dynamics of the gravitational system. 

We can estimate the length scale £ that characterizes 
the strength of thermophoretic interactions from the ex- 
periment of Ref. Q. For R = 1 /im, we can estimate 
that for a fully coated bead that is not self-propelled due 
to lack of asymmetry (Pe = 0), £ ^ 10 /im, while for 
the self-propelled colloids we could have a reduction by 
two orders of magnitude, namely, ^ ^ 0.1 /im. Consider- 
ing the confinement length to be /i ^ 10 — 100 /im, we 
find that it is very easy to realize a sufficiently dilute ex- 
perimental system, which is in the strong coupling limit. 
While the laser intensity provides a continuous tuning pa- 
rameter, the presence or absence of self-propulsion could 
move the system much faster in the parameter space. We 
note that Eq. ^ could be used to study the time depen- 
dence of the nonlinear dynamics of the colloids as in the 
analogous electrokinetic system [lit . 

The are a number of effects that we have not considered 
in the present analysis. We have neglected the temper- 
ature dependence of Deff in the calculation that led to 
Eq. ([6]), which could introduce corrections of the order 
of AT/T to the argument of the exponential. Moreover, 
hydrodynamic interactions have been shown to lead to 
nonlocal relations between the temperature profile and 
the diffusion coefficient of tracer particles [31 . However, 
we do not expect these effects to change the qualitative 
behavior of the system. We have also neglected the hy- 
drodynamic interaction between the colloidal particles 
themselves. This is justified, as sclf-thermophoretic col- 



loids are effectively source dipoles and lead to velocity 
fields that decay as 1/r^, which is faster than the ther- 
mophoretic interaction that decays as l/r^. 

In conclusion, we have studied the collective behavior 
of active colloids that act as mobile heat sources, and 
found that thermo-repulsive colloids could organize into 
hollow bands, tubes, or shells, depending on the geom- 
etry, while thermo-attractive colloids could go unstable. 
We note that similar equations could be used to study 
collective chemotaxis of diffusiophoretically active parti- 
cles using the analogy of nonequilibrium phoretic phe- 
nomena [7[. 
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